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PROCEDE ET SYSTEMS DE MESURE PES DEGRADATIONS D'UNE 
IMAGE VIDEO IN TRQDUITES PAR UN CODAGE A REDUCTION DE 
DEBIT. 

6 

La pr6sente invention conceme un proc^d6 et un systfeme de mesure des 
degradations d'une image vid^, inti-oduites par un systeme de codage k 
reduction de d6bit. 

1 0 Elle s'applique notamment, mais non exclusivement, au domaine d^ r^seaux de 
diffusion de signaux audiovisuels num^riques k has d6bit ou h trd's bas d^bit, 
ainsi qu'au domaine de la production de tels signaux. Elle s'applique en 
particulier k la surveillance de la quality de service d'un r6seau de diffusion de 
signaux audiovisuels num6riques. 

16 

. La num6risation des signaux video permet de copier, stocker et transmettre ce;. 
type d' informations en maintenant une quality d' image constante; Cependant, laf 
grande quantite d'informations v6hicul6e par les images vid^o n6cessite eii: 
pratique I'utilisation de m6thodes de compression num^rique pour r^uire lei 

20 d6bitbinaire. ^ 

'ft- 

Une methode de compression trhs r^andue dans le domaine de la vid6o est:; 
d^crite dans la norme ISO/CEI 13918 MPEG2. Get algorithme. est dit de typi 
"avec pertes" car I'image restitu6e aprds d^codage n'est pas identique k 
25 I'original. Get algorithme est has6 sur un d6coupage de I'image en blocs et sur 
r application d'une transform^e par exemple de type transform^ en cosinus 
discrete aux pixels de chaque bloc, ce qui permet d'obtenir une representation 
fr6quentielle de I'amplitude de la luminance des pixels sous la forme d'autant 
de coefficients que de pixels dans le bloc. 

30 

Afm de maintenii- une qualit6 acceptable pour le tei6spectateur final, les 
algorithmes de compression tiennent compte des propri^t^s de perception du 
systdme de vision humain. Cependant, les contraintes de debit impos^es par les 
systemes de transmission n6cessitent I'application de taux de compression qui 
35 influent sur la quality de I'image per9ue par le teiespectateur. 

II s'avdre que I'iraportance des degradations engendrees par le codage depend a 
la fois du taux de compression et de la complexite des images. Ces degradations 
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sont d*autant plus importantes que T image est complexe en ce qui conceme 
notamment le mouvement d'objets, la luminosite, et la texture. 

Parmi les degradations qui apparaissent dans les images h la suite de 
5 Tapplication de la technique de codage MPEG2, on peut citer les erreurs 
granulaires, les deformations de contoxirs, T apparition de contours dit 
"exotiques*' et les effets de blocs. 

II apparatt done n^cessaire d'6valuer de fa9on permanente la quality des images 
10 diffusees. II existe k cet effet des m6thodes d' Evaluation subjectives, largement 
utilis6es, faisant appel h TapprEciation humaine. Ces m^thodes sont toutefois 
lourdes h mettre en oeuvre et ne peuvent pas 6tre utilis6es en temps r6el sur un 
feseau de diffusion en fonctionnement^ 

II existe Egalement des m6thodes dites "avec reference" bas6es sur la 
comparaison de T image dent on veut 6valuer la qualite avec une image de 
reference. L'image de r6f6rence est en general celle qui correspond a Timage k 
analyser avant son codage et/ou sa transmission. Cette solution s'avere peu 
pratique du fait qu'elle necessite un acces a une ou plusieurs images de 
reference. En outre, dans le cas oil T image video est transmise, il se pose le 
probldme du transport de T image de r^f^rence vers le lieu de reception de 
rirhage k analyser. 

D'autres solutions dites "sans r6f6rence" permettent d' analyser 
25 automatiquement des images sans avoir a effectuer ime comparaison avec des 
images de r6f6rence. L*efficacit6 et la robustesse de chacune de ces solutions 
resident dans la methode employ6e pour mesurer des paramfetres relatifs k la 
quality de T image, 

Certaines de ces solutions sont bas6es sur une detection de T effet de bloc 

30 effectu6e dans le domaine spatial par des calculs de gradients aux frontidres des 

blocs. Pour 6viter de confondre la fronti&re d*un objet de Timage avec un effet 

de bloc, le gradient est compart avec des gradients intra bloc. L*effet de bloc est 

d6tect6 a Taide d*un crit^re de decision applique au comportement des gradients 

» 

inter et intra bloc. 

35 

Ainsi, la soci6t6 Rohde et Schwarz a d6velopp6 une m6thode de detection 
d* effet de blocs consistant k calculer pour chaque macro bloc de l'image un 
vecteur de gradient horizontal, et k calculer ime moyenne sur toute Timage de 
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chaque composante du vecteur. Les variations au cours du temps des 
composantes de ce vecteur permettent de mettre en 6vidence des composantes 
au comportement marginal qui repr&entent les frontidres de blocs d6grad6es 
par la transform^e du traitement de compression. La detection de ces 
5 composantes marginales permet de determiner un critfere de d^ection d'effet de. 
bloc repr^sentatif de la degradation de I'image, 

Ce principe de calcul de gradient est 6galement d6crit dans le brevet FR 2 805 
429 d6pos6 par la Demand^esse. Cette demande de brevet d^crit un proc^dd 
10 base sur la combinaison d'une image de gradient binaris^e et d'june image de 
"pseudo vecteurs" de mouvement, calcul6e k partir d'au moins deux images 
successives. La combinaison de ces deux images permet d'estimer un taux de 
faux contours dans I'image, utilise ensuite pour ^valuer une note de qualite. 

1 5 Dans le brevet FR 2 785 1 16 depose par la Demanderesse, les gradients calculus 
sur toute I'image a analyser sont soumis k des filtres psychovisuels 
traduisent I'effet de masquage contextuel. Un taux de fronti^res de blocj 
visibles detectes est ensuite calcul6 en recherchant une pseudo periodicite parmi 
les gradients de forte valeur, la quality de I'image etant evalu^e sur la base de ce 

20 taux. 

II s'avfere que les m6thodes basees sur le calcul de gradients appliquent des 
filtres pour ne s'interesser qu'^ un certain type de contenu de I'image : 
frontieres ou hautes frequences. Ces methodes ne permettent done d'analyser 

25 qu'une partie de I'information contenue dans I'image. II en resulte que leur 
fiabilite en terme de detection des degradations de I'image se trouve limitee. En 
outre, les precedes bases sur I'usage des gradients pour I'estimation de frontiere 
ou la detection de contours sur, I'image, sont relativement sensibles au bruit, ce 
qui afifecte la fiabilite de I'estimation de la qualite du contenu intrinseque de 

30 I'image. 

Par ailleurs, les methodes basees sur le calcul de moyenne sur I'image entidre 
reduisent de manidre drastique I'importance des degradations locaUsees dans 
une partie de I'image, ce qui rend difficile la detection de telles degradations 
35 localisees, et done affecte la fiabilite de revaluation de la qualite de I'image. 

Pour reduire ces inconvenients, certaines methodes prevoient d'analyser 
plusieurs images successives. Ces methodes ne sont done pas exploitables pour 




-4- 



analyser une image isol6e hors du cadre de la vid6o. 

La pr^sente invention a pour but de supprimer ces inconv6nients. Cet objectif 
est atteint par la provision d^in proc6d6 de mesure des degradations d'lme 
5 image numerisee, introduites par un codage de Timage consistant k decomposer 
r image en blocs de codage suivant une grille de codage et k appliquer sur des 
doim^es de pixels de chaque bloc un traitement de codage faisant intervenir un 
calcul de transformee par blocs et un calcul de transform^e par bloc inverse. 
Selon r invention, ce procede comprend des 6tapes consistant k : 

10 

- determiner la grille de codage de I'image cod6e, afm de retrouver la 
decomposition en blocs de codage de Timage, utilis6e lors du codage de 
I'image, 

- d6caler la grille de codage par rapport k Timage cod6e, de manifere a definir 
15 une decomposition de Timage en blocs d' analyse couvrant chacun une 

frontifere entre deux blocs de codage adjacents, 

- appliquer le calcul de transformee par bloc aux donnees de pixels de T image 
codec en utilisant la grille de codage decaiee pour obtenir des coefficients 
transformes pour chaque bloc d'analyse defini par la grille de codage 

20 decaiee, 

- extraire des coefficients transformes, des coefficients susceptibles d'etre 
affectes par un effet de bloc resultant du codage de T image, 

- appliquer aux coefficients transformes extraits le calcul de transformee par 
bloc inverse pour determiner des donnees de pixels de chaque bloc d' analyse, 

25 - estimer pour chaque bloc d' analyse un indicateur de degradation due aux 
effets de bloc, a partir des donnees de pixels de Timage codec et des donnees 
de pixels de chaque bloc d'analyse, obtenues par le calcul de transformee 
inverse, et 

- determiner une mesure de degradation de Timage en faisant la somme des 
30 indicateurs de degradation de chaque bloc d' analyse. 

Selon une particularite de Tinvention, I'estimation d'un indicateur de 
degradation pour chaque bloc d^analyse comprend les etapes consistant k : 

35 - calculer a partir des donnees de pixels obtenus pour le bloc d' analyse une 
moyenne de differences inter pixel au niveau de la frontiere inter bloc de la 
grille de codage, couverte par le bloc d' analyse, 

- calculer une moyenne et im ecart type portant sur les pixels des deux blocs 
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adjacents de la grille de codage, pardellement couverts par le bloc d'analyse, 

- calculer un facteur de pondaration en fonction de la moyenne et de I'^cart 
type obtenus pour le bloc d'analyse, et 

- calculer une activity spatiale du bloc d'analyse k partir d'activit^s spatiales 
daermin&s pour chacun des deux blocs adjacents de la grille de codage 
partiellement couverts par le bloc d'analyse, 

- I'indicateur de degradation du bloc d'analyse 6tant determine en fonction de 
la moyenne calcul6e des diff^nces inter pixel, du facteur de pond6ration et 
de I'activita spatiale du bloc. 



Avantageusement, I'indicateur de degradation d'un bloc d'analyse est obtenu k 
Taide de la formule suivante : 



dans laquelle Alij est la moyenne de differences inter pixel au niveau db4& 
frontidre inter bloc de la grille de codage, couverte par le bloc d'analyse, wy ^est 
le facteur de pond^ration, \|/ est une constante pred6finie, et ACTy est l'activit6 
spatiale du bloc d'analyse. -^^ 



Selon une autre particularite de 1' invention, le calcul de transformee est 
20 applique aux blocs de codage de I'image codee, les activites spatiales 
determinees pour chacun des deux blocs de codage etant obtenues a partir des 
coefficients transformes de chacun des deux blocs de codage. 

Avantageusement, les activites spatiales determinees pour chacun des deux 
25 blocs de codage sont obtenues k partir des formules suivante : 



dans lesquelles DCy et Acj;j(u,v), et respectivement Dcg et Ac3(u,v) sont 
les coefficients transformes de chacun des deux blocs de codage adjacents 
30 partieUement couverts par le bloc d'analyse, et NiU(u,v) est une fonction de 
masquage modeiisant un masquage par voisinage. 



10 





Acg<u,v)NiU(u,v) 



AcS(u,v)Nill(u,v) 
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De pr6f6rence, la moyenne et I'^cart type calculus pour chaque bloc d' analyse 
sont d6termin&5 partir des coefficients transformes de chacun des deux blocs 
de codage adjacents partiellement couverts par le bloc d' analyse. 

5 Selon une autre particularite de I'invention, le facteur de pond^ration est obtenu 
par la formule suivante : 



20 



25 



/ 



1+ 



si mjr^C 



Inl 1+-*: ^ I sinon* 



dans laquelle 



1+ 



'J255-C, 



l + CT; 



Infl+r^ 



Hij et ay sont respectivement la moyenne et I'^cart type calculus pour chaque 
10 bloc d'analyse, et C, est un paramdtre correspondant au maximum de la 
sensibility de Toeil humain. 

Selon encore une autre particularity de Tinvention, on effectue une s61ection des 
blocs d'analyse susceptibles de contenir un efifet de bloc, pr^alablement k 
1 5 r estimation d*un indicateur de degradation pour chaque bloc d'analyse. 

Avantageusement, la selection pr^alable comprend une 6tape consistant k 
6carter les blocs d'analyse dont les coefficients de transform6s extraits sont 
sup6rieurs k un seuil pr6d6fini. 



De pr6f6rence, la s61ection prealable comprend une etape consistant k 
s61ectionner les blocs d'analyse ayant, au niveau de la firontidre inter bloc de la 
grille de codage, couverte par le bloc d'analyse, des pixels pr6sentant une 
ynergie repr6sentant une proportion significative de 1' Anergic du bloc. 

Selon encore une autre particularity de Tinvention, le d^calage de la grille de 
codage est effectu6 horizontalement par rapport k I'image cod^e. 



Altemativement, le d^calage de la grille de codage est effectu^ verticalement 
30 par rapport k V image cod^e. 
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De pr6f6rence, le calcul de transfomi^e par blocs est un calcul de tfansform^e 
en cosinus disor^te. 

L'invention conceme 6galement un systdme de mesure des degradations d'une 
image mmi6ris6e, introduites par un codage de I'image, comprenant des moyens 
de calcul pour mettre en oeuvre le proc6d6 ddfini ci-avant. 

Un mode de realisation pr^f^r^ de l'invention sera d6crit ci-apr^s, & titre 
d'exemple non limitatif, avec r6f6raice aux dessins annex6s dans lesquels : 

La figure 1 repr&ente scMmatiquement le systdme de mesure 
selon l'invention, int6gr6 dans une chaine de traitement d'image ; 

La figure 2 montre plus en d6tail le syst^me de mesure repr6sent6 
sur la figure 1 ; 

La figure 3 montre deux blocs de pixels de I'image pour illustrer 
le proc6d6 selon l'invention ; a 

La figure 4 montre plus en detail une partie du syst&tne tepr^sent^*^' 
sur la figure 2 ; ^- 

La figure 5 montre une courbe repr&entant une fonction de 
pond^ration utilise par le ptoc6d6 selon l'invention ; 

La figure 6 montre d'autres courbes illustrant ime autre fonction 
de pond6ration utilis6e par le proc6d6 selon l'invention ; 

La figure 7 montre des courbes de variation de la mesure de 
degradation obtenue conform6ment au proc^de selon l'invention, 
appliqu6e k une sequence d' images. 

La figure 1 repr^sente un syst&ne de mesure 1 de degradations selon 
l'invention, qui est destine k 6tre utilise comme une sonde susceptible d'Sti-e 
appliquee en tout point d'une chaSne de traitement d'images numeriques, telle 
qu'une chaJne de diffusion ou de production de signal video. 




-8- 



Le systfeme de mesure 1 est plus particulidrement con9U pour mesurer la qualit6 
des images en sortie d'un systeme de codage 10 h reduction de ddbit respectant 
la norme MPEG 2. A cet effet, il se base sur une combinaison de filtrage 
fir^quentiel et temporel. 

5 

Ce systdme est applicable h chaque fois qu'il est n6cessaire d' identifier des 
d^fauts de codage d'un signal vid6o num6rique, notamment afin de determiner 
le d6bit convenable pour une sequence d' images donn6e, en fonction de la 
quality attendue. 

10 

Le principe de codage selon la norme MPEG 2 consiste k decomposer Timage 
num6ris6e en blocs de N x N pixels (N etant par exemple 6gal k 8) selon ime 
grille de codage, et d'appliquer h chaque bloc une transfonn6e en cosinus 
discrete (DCT : Discrete Cosine Transform) permettant de passer du domaine 
15 spatial au domaine fr6quentiel, puis d'annuler certaines composantes de la 
transform6e correspondant h des frequences eiev6es, avant d'appliquer une 
transformation DCT inverse pour repasser dans le dommne spatial, c'est-^-dire 
recuperer les pixels du bloc correspondant, 

20 Le traitement de transformee en cosinus discrete consiste a appliquer k chaque 
bloc Bij- de N X N pixels de I'image un calcul de coefficients de transformee AC 
et DC obtenus k I'aide des formules suivantes : 

DCi,j=F^j(0,0)=;^^'^fy(x,y) (1) 

et ACy(u,v)=Fg(u,v) , avec u + v 0 (2) 

25 dans lesquelles : 

fy(x,y) repr^sente la luminance du pixel au point (x, y) dans le bloc By, x et y 
etant les index de position respectivement horizontale et verticale du pixel dans 
le bloc Bij, et u et V sont compris entre 1 et N-1 et repr^sentent respectivement 
30 les index de la frequence spatiale horizontale et verticale, respectivement, et 

c(0) =-4= et c(u) = 1 si u 96 0. 
v2 

Le traitement inverse de la transformee en cosinus discrete consiste k appliquer 
aux coefficients DQ j et ACy(u,v) de chaque bloc Bjj la formule suivante : 




Sur la figure 2, le systdme de mesure 1 comprend un module de synchronisation 

11 permettant de synchroniser chaque image cod6e k trdter avec sa grille de 
5 codage, c'est-i-dire la d6composition de I'image en blocs de N x N pixels. Un 

exemple de traitement de synchronisation susceptible d'etre effectu6 par ce 
module est ddcrit dans la demande de brevet FR 2 769 452 d6pos6e par la 
Demanderesse. Ce traitement de synchronisation permet en fait de rep6rer dans 
I'image la decomposition en blocs utilise par le traitement precedent de codage 
10 effectue par le systdme 10, et done de determiner k quel bloc de codage By 
appartient chaque pixel de Pimage. 

L'image codec et sa decomposition en blocs sont appliques a deux branches de 
traitement 2, 3, la premidre branche 2 comprenant un module 16 de transform^ 
15 en cosinus discrete, et la seconde branche 3 comprenant tout d'abord un module 

12 de decalage horizontal d'un demi bloc de la grille de codage par rapport -4 
I'image. La figure 3 illustre ce decalage horizontal d'un demi bloc, en montrant 
un bloc By de 8 X 8 pixels d^caie horizontalement d'un demi bloc par rapport a 
la grille de codage illustree par les blocs 2 1 et 22 adjacents, 

20 r 

Dans la seconde branche de traitement 3, I'image codec et la grille de decodage 
decaiee en sortie du module de decalage 12 sont appliquees h I'entree d'un 
module 13 de transformee en cosinus discrete. Un module d'extraction 14 de 
frequences horizontales traite ensuite les coefficients de transformee issus du 
25 module 13. Ce module consiste en fait k ne conserver dans le domaine 
frequentiel que. les firequences (coefficients ACij(u,v) et DCy) susceptibles de 
rendre visible un effet de bloc vertical. Dans le cas ou N=8, on peut ainsi 
estimer que doivent dtre conserves uniquement les coefficients DCy, ACij(u,v) 
tel que ; 

30 0 ^ u < 8 si V = 0, 

2^u<8 siv= 1, et 
6^u<8siv = 2. 

Les coefficients conserves ACij(u,v) et DC,- j sont ensuite traites par un module 
35 15 de transformation DCT inverse pour recuperer les informations spatiales 
I|j(x,y) = I(8i+x, 8j+y) relatives aux frequences conservees, obtenues k I'aide de 
la formule (4). On peut k cet effet se limiter k la quatrieme et la cinquieme 
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colonne 24, 25 concem^es par I'effet de bloc du au codage, et qui 
correspondent k la j&ontidre 23 entre les blocs 21 et 22 de la grille de codage sur 
ia figure 3. 

5 En sortie de la branche 3, on dispose done des domi6es de pixel I(8i+x, 8j+3) et 
I(8i+x, 8j+4) de toute I' image, i et j etant les index des blocs decal6s de P image 
selon la grille de codage d6cal6e obtenue en sortie du module de d^calage 12. 
En sortie de la branche 2, on dispose des coefficients DCy et ACy(u,v) de 
chaque bloc de I'image relatifs h la grille de codage en sortie du module de 

10 synchronisation 11. Pour chaque bloc By d6cal6, on dispose des coefficients 
DCyetACy(u,v) du bloc de gauche et 21 des coefficients DCy et 

Ac2(u,v) du bloc de droite 22. 

Les donn6es issues des deux branches de traitement 2 et 3, c'est-a-dire en sortie 
15 des modules 16 et 15, sont appliqu6es en entr6e d'un module 17 d'6valuation de 
I'effet de bloc. 

Sur la figure 4, les donn^es en entree du module 17 sont appliquees aux 
modules suivants : 



- un module 31 de calcul de moyenne des differences inter pixels k la fronti^>re 
verticale 23 de chaque paire de blocs adjacents 21, 22 de I'image cod6e ; 

- un module 32 de calcul de moyenne et d'^cart type de la luminance de 
diaque bloc, et 

25 - un module 34 de calcul de contraste local. 

Pour 6valuer la moyenne des diff6raices inter pixels sur chaque bloc d€cal6 By, 
le module 31 applique aux donn6es de pixel issues de la seconde branche 3 la 
formule suivante : 



La moyenne ainsi calcul^e d€crit en fait le comportement de l'op6rateur 
gradient k la fronti^re verticale 23 de deux blocs adjacents 21, 22 issus de la 
grille de codage de I'image. 




k=0 



(5) 



35 Dans le module 32, une moyenne Hy et un 6cart-type ay de luminance sont 
estimes pour chaque bloc By sur les deux blocs adjacents 21, 22 correspondants 




de la grille de codage. On pourra par exemple utiliser les coefficients AC et DC 
des blocs de codage dans le domaine DCT foumis par la premiere branche 2 
selon les formules suivantes : 



Kj=i(DCa+DcPj) 



et 



V u+v?iK) ■ u+v^ 



(6>, 



(7) 



Bien entendu, la moyenne et I'^cart type ay peuvent dgalement atre obtenus 
1 0 directement k partir des pixels des blocs 21 , 22 correspondant a chaque bloc By, 
k I'aide de formules classiques de calcul de moyenne et d'feart ^e. 



15 



Les moyennes et hearts types calculus par le module 32 sont ensuite appliques t 
un module 33 de calcul d»un facteur de pond^ration Wy(ny, cyy, Q pour chaque 
bloc By. Pour calculer ce facteur de pond^ration, le module 33 applique la 
fonction de pond^ration d€finie par la fonnule suivante : 



20 



25 



1+ 



hi 



si Hy^^ 



In 



Gil 



(8j' 



smon> 



1+cr 



(9) 



et C, est un paramdtre correspondant au maximum de la sensibility de I'oeil. 
Cette fonction de pond^ration est d6crite dans la publication "A Generalized 
Block-Edge Impairment Metric for Video Coding" de H.R. Wu et M. Yuen, 
ffiEE Signal Processing Letters, Vol. 4, N^ll, Novembre 1997. Toutefois, la 
fonction de ponddration d6crite dans cette publication pr6sente une 
discontinuity au point = C lorsque I'^cart type ay est non nul. Cette fonction 
a done €t6 adapt^e par le proc^dg selon I'invention pour atre continue en ce 
point. 
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La figure 5 illustre cette fonction de ponderation sous la forme d'une courbe 41 
en fonction du niveau de gris ou de la luminance moyenne fiy des pixels du 
bloc, dans le cas ou ^ est choisi 6gal a 81, ce qui correspond au maximum de la 
5 sensibilite de I'oeil humain, et aix o est 6gal h 20. 

Le calcul du contraste local effectu6 par le module 34 consiste h estimer 
ractivit6 spatiale ACTy de chaque bloc By dans le domaine transform^, h 
partir des coefficients DCij et ACy, les coefficients ACy €taaA pond6r6s par une 
10 fonction de masquage pour mod^liser le masquage par le voisinage. L' activity 
spatiale ACTy de chaque bloc By est ainsi estim^e en appliquant la formule 
suivante : 

ACTy =^(ACt3 + ACli^) (10) 
oil ACtS= i; rACy(u,v)Nill(u,v)l (11) 

ACT-^— ^-^J t rAc5(u,v)Nill(u,v)f , (12) 
l+DCyVu.v=oL " -I 

U+V940 

dans laquelle Nill(u,v) est la fonction de masquage. qui est par exemple celle 
decrite dans le document "A Visual Model Weighted Cosine Transform for 
Image Compression and Quality Assessment", de N. B. Kill, IEEE Transactions 
on Communications, Vol. COM-33, N° 6, Juin 1985. 



15 et 



20 



Dans ce document, la fonction de masquage, est donn6e par la formule 
suivante : 

Nill(u,v)=A(©)H(co) (13) 

ou oi)=Vu^+v^ represente la frequence radiale exprim6e en cycles par degr6 (u 
25 etv e[0,N-l]), 



A(a))=J|+4 



2 

n 



( r~2 

1 « V a 



2 

f-+l 
-0.290) 



(14) 



H(a))=(0.31+0.69co)e (15) 
et a est commun6ment 6gal k 1 1,636 deg'' pour un traitement de la luminance 
sur 8 bits. 



.30 
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La figure 6 illustre cette fonction de pond6ration en repr6sentant en fonction de 
la frequence spatiale ©, la courbe 43 de variation de H((o), la courbe 44 de 
variation de A(<d) et la courbe 45 de variation de la fonction Nill(a>). 

5 Dans les formules (1 1) et (12), la normaUsation par le coefficient DC permet de 
consid^rer Fad^quation au contraste local, conform6ment k la loi de Weber, et 
I'unit^ est ajout^e au d6nominateur pour que la fonction soit continue en z6ro. 

Les r^sultats des calculs effectu6s par les modules 31, 33 et 34 sont enfin 
10 appliques k un module 35 d' Evaluation de la degradation o de I'image du fait de 
Teffet de bord resultant du codage. Pour estimer une mesure de, degradation 
d'image, ce module effectue le calcul suivant : 

"=Ii«y (16) 



15 repr^sente Timpact des efifets de bloc sur le bloc By et x^t est une constanfe 
daermin^e expErimentalement qui permet d»ajuster Pimportance attribute 'k 
I'activite spatiale par rapport k la moyenne pond6r6e des differences inter 
pixels. ^ 

En fait, la mesure de degradation o estim6e par ce module correspond k un 
20 indicateur de visibilite par Toeil humain des degradations dues k I'effet de bloc- 
La figure 7 montre sous la forme de courbes 46, 47, 48, les variations de cette 
mesure en fonction du numero d'image dans une sequence d' images video, 
respectivement pour des debits de 9, 6 et 4 Mbits/s. 

25 

Ces courbes montrent que les degradations des images sont amplifiees au fur et 
k mesure que I'on reduit le debit. Le codage grossier par quantification des 
coefficients DCT provoque done Papparition de I'effet de bloc. La partie de la 
sequence d'images oil la mesure de degradation est eievee correspond k des 
30 images dont le contenu spatial et temporel evolue brutalement. 

Pour optimiser les traitements effectues par le module 17, on pent soumettre les 
donnees de bloc en entree k une selection k I'aide d'un module de selection des 
blocs Bij susceptibles de contenir un effet de bloc. A cet effet, ce module 
35 compare I'energie des colonnes centrales 24 et 25 de chaque bloc By avec 
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r6nergie du bloc. Le bloc est retenu pour restimatlon de la d6gradation si 
r6nergie des coloimes 24 et 25 repr^sente line proportion significative, par 
exemple sup^rieure d 25%, de r6nergie du bloc. Avant d'effectuer cette 
selection, on peut 6galement ^carter les blocs dont les coefficients DCT retenus 
5 par le module d' extraction 14 sont inf^rieurs a un seuil pr6d6fini, par exemple 
voisin de z6ro. 

II est k noter que le traitement qui vient d'etre dterit est con9u pour d^tecter un 
effet de bloc vertical. Bien entendu, ce traitement peut 6tre modifi6 d'une 

10 maniere 6vidente pour d6tecter un effet de bloc horizontal. II suffit pour cela 
d'effectuer dans le module 12 un d^calage vertical d'un demi bloc, puis de 
s61ectionner dans le module 14 d'extraction de frequences, les frequences 
suscqptibles d*6tre affectees par un effet de bloc horizontal, et enfin de ne 
conserver parmi les pixels reconstitufe dans le module 15 que les pixels des 

15 quatriemes et cinquidmes lignes du bloc d6cal6 correspondant k la frontifere 
horizontale entre les deux blocs non d6csl6s adjacents, correspondant au bloc 
d6cal6. 

On pourrait 6galement effectuer en parall^le des traitements de detection 
20 d'effets de bloc verticaux et horizontaux et d'utiliser conjointement les r^sultats 
foumis par ces deux traitements dans le module 35 par exemple. 

Apr^s comparaison des mesures objectives obtenues par le proc^de et le 
systeme selon I'invention. aux notations effectu6es lors de tests subjectifs par 
25 des panels de personnes representatives, on peut remarquer une forte correlation 
entre les mesures objectives et les mesures subjectives, ce qui assure de 
I'efScacite du procede et systeme selon I'invention. 




REVENDICATIONS 

1. Proc6d6 de mesure des degradations d'une image num6ris6e, 
introduites par un codage de I'image consistant h decomposer T image en .blocs 

5 de codage suivant une grille de codage et k appliquer sur des donn^es de pixels 
de chaque bloc im traitement de codage faisant intervenir un calcul de 
transfonn6e par blocs et un calcul de transformee par bloc inverse, 
caract^rise en ce qu'il comprend des 6tapes consistant k : 

- determiner (11) la grille de codage de 1' image codde, afin de retrouver la 
10 decomposition en blocs de codage de Timage, utilisee lors du codage de 

rimage, 

- decaler (12) la grille de codage par rapport a I'image codee, de manidre h 
definir une decomposition de I'image en blocs d'analyse (Bij), couvrant 
chacun une frontidre (23) entre deux blocs de codage (2 1 , 22) adjacents, 

15 - appliquer (13) le calcul de transformee par bloc aux donnees (fy(x,y)).de 
pixels de Timage codee en utilisant la grille de codage decaiee pour obte^ 
des coefficients transformes (DCy, ACy(u,v)) pour chaque bloc d'analyse 
(By) defini par la grille de codage decaiee, 

- extraire (14) des coefficients transformes, des coefficients susceptibles d'6jre 
20 affectes par uneffetde bloc resultant du codage de rimage, .iJ^ 

- appliquer (15) aux coefficients transformes extraits le calcul de transformee 
par bloc inverse pour determiner des donnees (Iij(x,y)) de pbcels de chaque 
bloc d'analyse, 

- estuner (17) pour chaque bloc d'analyse un indicateur (uy) de degradation 
25 due aux effets de bloc, partir des donnees (fy(x,y)) de pbcels de I'image 

codee et des donnees de pbcels aij(x,y)) de chaque bloc d'analyse. obtenues 
par le calcul de transformee inverse, et 

- determmer (17) une mesure de degradation (u) de I'miage h partir des 
indicateurs de degradation (uy) de chaque bloc d'analyse. 

30 

2. Procede selon la revendication 1, 

caracterise en ce que I'estimation d'un indicateur de degradation (oy) pour 
chaque bloc d'analyse (By) comprend les etapes consistant k : 

- calculer k partir dea donnees (Iy(x,y)) de pbcels obtenus pour le bloc 
35 d'analyse une moyenne de differences inter pixel (Aly) au niveau de la 

frontiere (23) inter bloc de la grille de codage, couverte par le bloc d'analyse, 

- calculer une moyenne (fxy) et un ecart type (ay) portant sur les pixels des 
deux blocs (21, 22) adjacents de la grille de codage, partiellement couverts 
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par le bloc d' analyse, 

- calculer un facteur de pond6ration (Wij) en fonction de la moyenne et de 
r€cart type obtenus pour le bloc d' analyse, et 

- calculer une activity spatiale (ACTy) du bloc d'analyse a partir d'activit6s 
spatiales (ACTy , ACTy) determin6es pour chacun des deux blocs (21, 22) 
adjacents de la grille de codage partiellement couverts par le bloc d' analyse, 

Tindicateur de d6gradation du bloc d*analyse 6tant d6termin6 en fonction de la 
moyenne calcul6e des differences inter pixel, du facteur de ponderation et de 
ractivit6 spatiale du bloc. 
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3. Proc^d6 selon la revendication 2, 
caracteris6 en ce que Pindicateur de degradation d'un bloc d' analyse est obtenu 
a I'aide de la fonnule suivante : 

'•J l+V,/|ACTy| 

15 dans laquelle Alij est la moyenne de differences inter pixel au niveau de la 
firontiere (23) inter bloc de la grille de codage, couverte par le bloc d'analyse, 
Wy est le facteur de pond6ration, \\i est une constante predefinie, et ACTy est 
r activity spatiale du bloc d'analyse. 

20 4. Procede selon la rev^dication 2 ou 3, 

caracterise en ce que le calcul de transformde est applique aux blocs de codage 
de r image codec, les activites spatiales (AClfj,ACTy) determinees pour 
chacun des deux blocs de codage (21, 22) etant obtenues h partir des 
coefficients transformes (DCy, ACy(u,v)) de chacun des deux blocs de codage. 
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5. Procede selon la revendication 4, 
caracterise en ce que les activites spatiales (ACTi j,ACTy) determinees pour 
chacim des deux blocs de codage (21, 22) sont obtenues h partir des formules 
suivantes : 



30 ACTi 

u+v^tO 



i; rAC^u,v)NiU(u.v)T , 




dans lesquelles DCj^j et ACj j(u,v) et respectivement DC^ et Acjj(u,v) sont 
les coefficients transfonn6s de chacun des deux blocs de codage adjacents (21, 
22) partiellement couverts par le bloc d'analyse, etNill(u,v) est une fbnction de 
masquage mod61isant un masquage par voisinage. • 

5 

6. Proced6 selon Tune des revendications 2 k 5, 
caract6ris6 en ce que la moyenne (jaij) et I'^cart type (oij) calculus pour chaque 
bloc d'analyse (By) sont daennin^s h partir des coefficients transfomi6s (DCij, 
ACij(u,v)) de chacun des deux blocs de codage (21, 22) adjacents partiellement 
1 0 couverts par le bloc d*analyse. 
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7. Proc&i6 selon I'une des revendications 2^6, 
caract6ris6 en ce que le fecteur de pond^ration (wy) est obtenu par la fonnule 
suivante : 

X\n 1+ 



oi^ 



V 

f 

In 1+ 
f 

In 14 



l+as; 

V255-C 



SI 



smon^ 



1+ai 



In 



1+ 



A. 



^ij et ay sont respectivement la moyenne et r6cart type calcules pour chaque 
bloc d'analyse (By), et C, est un param^tre correspondant au maximum de la 
sensibility de Toeil humain. 

8. Proc^6 selon Tune des revendications I hi, 

caract6ris6 en ce que pr^alablement a I'estimation d'un indicateur de 
degradation oy pour chaque bloc d'analyse (By), on efifectue une selection des 
blocs d'analyse susceptibles de contenir un effet de bloc. 

9. Proc6d6 selon la revendication 8, 

caract6ris6 en ce que la selection pr6alable comprend une 6tape consistant k 
^carter les blocs d'analyse (By) dont les coefficients de transform6s extraits 
sont superieurs k un seuil pr6d6fini. 
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10. Proced6 selon la revendication 8 ou 9, 

caract^rise en ce que la s61ection pr6alable comprend une 6tape consistant h 
s61ectionner les blocs d' analyse (By) ayant, au niveau de la firontiere inter bloc 
(23) de la grille de codage, couverte par le bloc d' analyse, des pixels (Iij(x,y)) 
5 pr6sentant une Anergic repr&entant une proportion significative de r6nergie du 
bloc. 

1 1 . Proced6 selon Tune des revendications 1 & 10, 

caract6ris6 en ce que le d^calage de la grille de codage est effectue 
1 0 horizontalement par rapport h Timage cod6e. 

12. Proc6d6 selon I'une des revendications 1 1 1, 

caract&is6 en ce que le d^age de la grille de codage est effectu6 
verticaleraent par rapport k Timage cod6e. 

15 

13. Proc6d6 selon Tune des revendications 1 a 12, 

caract6ris6 en ce que le calcul de transform^e par blocs est un calcul de 
transform6e en cosinus discrete. 

20 14. Systeme de mesure des degradations d'une image nuraeris6e, 

introduites par un codage de 1' image, 

caract6ris6 en ce qu'il comprend des moyens de calcul pour mettre en ceuvre le 
proc6d6 selon Tune des revendications 1 4 13. 
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